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Abstract 

In a previous study, the flavor-changing fermion-graviton interactions have been analyzed in 
the framework of the standard model, where analytical results for the relevant form factors were 
obtained at the leading order in the external fermion masses. These interactions arise at one-loop 
level by the charged electroweak corrections to the fermion-graviton vertex, when the off-diagonal 
flavor transitions in the corresponding charged weak currents are taken into account. Due to the 
conservation of the energy-momentum tensor, the corresponding form factors turn out to be finite 
and gauge invariant when external fermions are on-shell. Here we extend this previous analysis by 
including the exact dependence on the external fermion masses. Complete analytical results are 
provided for all the relevant form factors to the flavor-changing fermion-graviton transitions. 
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1 Introduction 



In a previous analysis pQ , following the study of Ref. [2j , we have discussed the structure of the per- 
turbative corrections to the graviton-fermion-antifermion (Tff) vertex in the Standard Model (SM), 
focusing our attention on the flavor diagonal sector. On the other hand, in [2] the one-loop electroweak 
corrections which generate the off-diagonal graviton-fermion-antifermion vertex, were computed at the 
leading order in the external fermion masses. These studies address the structure of the interactions 
between the fermions of the Standard Model and gravity, beyond leading order in the weak coupling, 
which have never been presented before in their exact expressions. The choice of an external (classi- 
cal) gravitational background allows to simplify the treatment of such interactions where the coupling 
is obtained by the insertion of the symmetric and improved energy-momentum tensor (EMT) into 
ordinary correlators of the Standard Model. 

We have addressed some of the main features of the perturbative structure of these corrections, 
presenting their explicit form, parameterized in terms of a certain set of form factors. We have also 
discussed some of their radiative properties with regard to their infrared finiteness and renormaliz- 
ability, the latter being inherited directly from the Standard Model, when the coupling of the Higgs 
to the gravitational background is conformal. 

In general, one expects that such corrections are small, although they could become more sizeable 
in theories with a low gravity scale [21 IU El El [8] . In particular, one can consider the possibility of 
including, in these constructions, backgrounds which are of dilaton type, with dilaton fields produced 
by metric compactifications. The same vertices characterize the interaction of a dilaton of a sponta- 
neously broken dilatation symmetry with the ordinary fields of the Standard Model QUI [HI 02] • 
This second possibility is particularly interesting, in view of the recent discovery of a Higgs-like scalar 
at the LHC. 

Perturbative studies of these vertices have their specific difficulties due to the proliferation of form 
factors, and the results have to be secured by consistency checks using some relevant Ward identities. 
These need to be derived from scratch using the full Lagrangian of the Standard Model, as discussed 
in [13] and pQ. In this study we are going to reconsider the gravitational form factor of a Standard 
Model fermion in the presence of a background graviton in the off-diagonal flavor case, which has been 
discussed before [2J, extending that analysis. One of the goals of this re-analysis is to include all the 
mass corrections to the related form factors, which has not been given before. These corrections are 
important in order to proceed in a follow-up work with a systematic phenomenological study. In this 
respect, mass corrections are important in order to extract the exact behavior of these form factors 
in the infrared and ultraviolet limits, which may be of experimental interest. We have compared our 
new results against the previous ones given in [2] in the limit of massless external fermions and found 
complete agreement. 

Our work is organized as follows. In section 2 we give the theoretical framework of the Standard 
Model Lagrangian in a curved space-time, assuming as a background metric the usual 4-dimensional 
one. In section 3 we discuss the technical details for the calculation of the flavor-changing fermion 
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matrix elements of the energy momentum tensor at the leading order in perturbation theory. The con- 
tributions coming from the counterterms of the wave-function renormalization is separately discussed 
in section 4. In section 5 we analyse the role played by the Ward identity from the conservation of the 
EMT, while in section 6 we give the results for the complete set of the relevant form factors entering 
in the flavor-changing matrix elements of EMT. Finally, in section 7 we give our conclusions. 



2 The Lagrangian 

We follow closely the layout of our previous work [1] where more details concerning the general 
structure of the action describing the coupling of the Standard Model to gravity can be found. We 
just recall, in order to make our treatment self-contained, that the interaction of the Standard Model 
fields with gravity is described by the action integral 

S = ~~2 j dAx "/~9 R + J d A x^f^gC S M (1) 
together with a term of improvement 

Si = xJ d A x\f^gR H^H (2) 

where R is the scalar curvature and H is the Higgs doublet. The identification of this second term goes 
back to |14| . The coefficient x is an arbitrary parameter which at the special value Xc = 1/6 renders the 
Lagrangian conformally symmetric when the scalar is massless and guarantees its renormalizability at 
the leading order in the gravitational k, where k 2 = WttGn, with Gn being the gravitational Newton's 
constant. For instance, in the case of the Higgs field, this takes place if we drop the quadratic terms 
in the Higgs potential. As in our previous work, our results are given for an arbitrary \- 

The Standard Model action Ssm is obtained by promoting the ordinary SM Lagrangian to a curved 
background, which is parametrized by the metric expansion g^ v = rj^ v + Kh^ u where rj^ u = (+, — , — , — ) 
and h^ u denotes the fluctuation of the graviton field around the flat limit. At this order the graviton- 
matter interactions, which we are going to evaluate in the flavor changing fermion sector, are described 
by Green's functions with a single insertion of the energy-momentum tensor 

(3) 

The complete Standard Model energy-momentum tensor includes several contributions which can be 
found in [T3] , 
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Ssm + Si 



3 The perturbative expansion 

The interaction of one graviton with two fermions of different flavor is summarized by the vertex 
function 

t^^zi^p^mpjJi) (4) 
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Figure 1: Diagrams of one-loop SM corrections to the flavor-changing graviton fermion vertex, where 
fij and pi j specify the fermion flavors and corresponding momenta respectively, with fi ^ fj. 

that we intend to study. Here pj (fj) and pi (fi) indicate the momenta (flavor) of initial and final 
fermions respectively. We will restrict to the case of flavor-changing transitions, namely fi ^ fj. In 
order to simplify the results we will also use the combinations of momenta p = Pi +Pj and q = Pj — Pi- 
The external states are taken on their mass shell, pf = mf and p'j = and can be either leptons or 
quarks. From now on, we will assume that nii ^ nij. In the last case, since the EMT is diagonal in 
color space, the color structure is rather trivial and therefore we omit it. 

At tree level the flavor-changing gravitational interaction is absent so that the leading order con- 
tribution comes from the quantum corrections. At one loop level, instead, we decompose the T^ v 
matrix element as 



where the first term on the r.h.s represents the pure vertex corrections induced by the gauge 
boson and its Goldstone (j)^ exchanges, while the last term, Tfy^, includes the usual counterterms 
(CT) coming from the wave-function renormalization insertions on the external legs. The inclusion 
of this last term Tq t is needed in order to get finite results for the matrix element T^ v , as it will 
be extensively discussed in section [H The finiteness of the result is just a consequence of the non- 
renormalization theorem of conserved currents, when applied to the case of a conserved EMT. 

We choose to work in the gauge where every massive gauge field is always accompanied by its 
unphysical longitudinal part. The diagrammatic expansion of is depicted in Fig IT] and is made of 
one contribution of triangle topology plus contact terms (see Fig. [1] (c) and (d)) with a fermion and 
a graviton pinched on the same external point. The Feynman rules are listed in Appendix [A] The 
computation of these diagrams has been performed in dimensional regularization using the on-shell 
renormalization scheme. To check the correctness of our results the Ward identity of the conservation 
of the EMT, which will presented in section has been verified explicitly. 

where fi j and pij specify the fermion flavors and corresponding momenta respectively, with fi ^ fj 
Due to the chiral V-A nature of the W interactions, we expand the flavor changing matrix element 
in terms of invariant amplitudes fk and tensor operators Ok as 




(5) 




(6) 
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with the tensor operators given by 



Of = (yy + yy ) P L Of = r?^ M_ 

Of = ( 7 /y + Y<f) Pl Of = p*y M_ 

Of = t?^ M+ Of = q^g" M_ 

Of = pfy M + Of = (p»q u + q^p u ) M_ ^ 

Of = g*y M+ Of = ^ ( 7 *y + yy) P R 



w 



Of = (p<y + g/y ) M + Of = ^ ( 7 *y + 7 «y ) P K 

where Pl,r = (1 T 7s)/2 and M± = rrijP R ± mjPx, and Uij(pij) are the corresponding fermion bi- 
spinor amplitudes in momentum space. 

This is the most general rank-2 tensor basis that can be built out of two momenta, p and g, a metric 
tensor and Dirac matrices 7^ and 7 s . Its expression has been given in [2]. 
For the form factors appearing in Eq.([6]) we use the following notation 

fk(p,q) = ^2>>f Fk(p,q,m f ) , (8) 
/ 

where we have factorized the term Xf = VfiVfj (the external fermions are assumed here to be quarks 
of down type), with Vij the corresponding CKM matrix element. 



4 Contribution from the wave- functions renormalization 

The matrix element corresponding to the vertex corrections is ultraviolet divergent and, due to the 
non-renormalization theorem of the conserved EMT, it is made finite by adding the contributions from 
the wave- function renormalization on the external legs, namely T<f • This last contribution can be 
easily determined by using the following method, as illustrated in pp. We promote the counterterm SM 
Lagrangian to a curved background and then extract in the usual way the appropriate renormalized 
Feynman rules for single insertions of the EMT on the fields of the Standard Model. The metric is 
taken to be flat after all the functional differentiations. Then, for the off-diagonal flavor contributions 
(i 7^ j) to T^, we have 



1 CT 



i^, w&mpjjj) = -uMl (yy + yy) {c^p L + c^p r 



+ 2 77^ 



C£- {m % P L - mj P R ) + Cg- (miP R - mj P L ) 



1.) 



Uj{pj) 



0) 



where 
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L± 



6Z%±5Z$ 



C 



R± 



2 K 6Z ij ±5Z v 



(10) 



with 5Z- j being the fermion wave function renormalization constants. In the on-shell renormalization 



scheme, which we have chosen for our computation, the renormalization conditions are fixed in terms 
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of the physical parameters of the Standard Model to all orders in the perturbative expansion. In 
particular for the fermion wave function renormalization constants with i ^ j one obtains 

2 



SZ 



mf — ttta 



Re< TYij ^^(rrij) + irii mj H^irrij) + (m i + nij) H^im^) 



SZ, 



mf — m- 

1 j 



Rel mj Sg(m^) + m, mj ZfAmj) + 2 m 3 - Sf (mj) ^ . 



(11) 



The symbol Re gives the real part of the scalar integrals in the self-energies but it has no effect on 
the CKM matrix elements. Its presence yields SZjj = bZ^ {m^ o m|J . Remember also that if the 

mixing matrix is real Re can obviously be replaced with Re. 

For completeness we give the Standard Model flavor changing self-energies (i ^ j) 



f 



mj + 2m 2 v ) Bi{p 2 , m^mw) + ™w 



where 



Z«(P 2 ) 



Blip 2 , ml, mf] 



G f 

• 47r a^2 m mj ^ VifV f j Bl (P 2 > m }> m w) 
E V *f V h m ) MP^rn), m 2 w ) , 



2 2 

m l ~ m o 



(12) 



(13) 



We have explicitly checked that the counterterm in Eq.Q is indeed sufficient to remove all the ultra- 
violet divergences of the matrix element so that is finite, as expected. 



5 The Ward identity from the conservation of the EMT 

The conservation of the energy-momentum tensor constraints the T^ u matrix element reducing the 12 
form factors defined above to a smaller subset of 6 independent contributions. We can derive the Ward 
identity by imposing the invariance of the 1-particle irreducible generating functional - which depends 
on the external gravitational metric - under a diffeomorphism transformation and then functional 
differentiating with respect to the fermion fields. We omit the details of this procedure, which has 
been discussed extensively in [13] and [1] for the TVV' and the Tff vertices respectively. The analysis, 
in this new case, follows similar steps. In momentum space, for the unrenormalized matrix element 
we obtain the Ward identity 

qpTft = UiipMp^ipJ-pVTijipj) + | (I^K* - aTyfo)) \ujipj) , (14) 

where a^ v = [7 M ,7^]/4 and Tij(p) is the fermion two-point function which is given by 



r ij(p) 



^ii{p 2 )i>PL + m P 2 )i>P R + SgCp 2 ) {rrn P L + m 3 P R ) 



(15) 
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The off-diagonal (in flavor space) two-point form factors T> L,R ' S (p 2 ) are explicitly given in Eq. (|12p . 
The renormalized Ward identity is instead much simpler than Eq. (|14p being just q^T^ = 0. It implies 
a set of homogeneous equations [2] for the renormalized form factors fk(p, q) 

p-qfi{p,q) + q 2 f2(p,q) = o, 
Mp,q) + q 2 f5(p,q)+p-qf6(p,q) + 7^rfi2(p,q) = o, 



p-qh(p,q) + q 2 f6(p,q) + 7rr2~f u ^P'^ = °. 



Mr 
• < 

p2 _|_ r ,2 



hip,q) + f7{p,q) + q fa(p,q) +p- qfio(p,q) - A 2 fnip^q) = o. 



h l />• <{) + !>■ 'i J si i>. <i ) + gVioCp, g) - 4m 2 fn(p,q) = 0. 



P-qfn(p,q) + q 2 fi2(p,q) = 0, (16) 

which provide a strong test on the correctness of our results and allow to reduce the number of 
independent contributions to the T^ v matrix element. 

6 Flavor-changing form factors 

In this section we present the explicit expressions of the renormalized form factors defined above. 
They have been computed in the on-shell case retaining the full dependence on the internal (rrif,mw) 
and external masses (rni,m,j) and on the virtuality, q 2 , of the graviton line. They are expressed in 
terms of the dimensionless ratios xs = (m 2 +m 2 )/q 2 , xd = (m 2 — m 2 )/q 2 , Xf = m 2 /q 2 , xw = fn\v/q 2 
and of the combination A = x D — 2xs + 1. We recall that mj is the mass of the fermion of flavor / 
running in the loop. 

Due to their complexity we expand our results onto a basis of massive one-, two- and three-point 
scalar integrals as 

7 

Fk(p,q,m f ) = Y / C l k I l (17) 



1=0 



where 



1*0 = 1, h= B (q 2 ,m 2 ,m 2 f ) , 

h = A (m 2 ) - Aoirrv^) , I 5 = B (q 2 , m 2 w , m 2 w ) , 

h = S (m|, mj, ra^-) , h = Co(m 2 , q 2 , m 2 , mj , , m 2 ^) , 



(18) 



I 3 = B (m 2 ,m 2 ,m 2 w ), I 7 =C (m 2 ,q 2 ,m 2 ,m 2 w ,m 2 ,m 2 ). 

We give the explicit results for the renormalized form factors F%, F3, F4, Fj, Fg and F\\ while the 
remaining six can be obtained exploiting the Ward identities derived in the previous section 

F2 = Xd Fx , 
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^5 = -^Fs + xlFt + ^-Fn, 
F e = -x D F 4 - r-^-i^ii , 

X D 1 2 X S X D 

Fg = 2—z-Fi ^F 7 + x D F s ^— Fn , 



F w = --^Fi-xoFs + ^-Fu 
Fn = -x D F u . (19) 



The coefficients Cl denning F^ in Eq. (|17p are given in the Appendix O Finally we remark that the 
i*3, F5, Fj and Fg form factors depend also from the parameter x which appears in the gravitational 
coupling of the (j)^ Goldstone bosons through the improved energy-momentum tensor. 



7 Conclusions 

We have presented the computation of the structure of the gravitational form factors of the Standard 
Model fermions in the off-diagonal flavor sector. The analysis has been developed according to our 
previous study pQ where we have discussed the electroweak corrections in the flavor conserving case. 
The work extends a previous investigation [2] of the same flavor-changing vertex in which the external 
mass dependence has not been included. The exact expressions presented in our work are relevant for 
a phenomenological study of the small and intermediate momentum behavior of these form factors, 
which we plan to address in the near future. 
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A Feynman rules 

We collect here all the Feynman rules involving a graviton that have been used in this work. All the 
momenta are incoming 

• graviton - gauge boson - gauge boson vertex 

ya 




(20) 

where V stands for the vector gauge boson W. 
graviton - fermion - fermion vertex 




*f (*i + W + l v (h + k 2 y - 2 rT (ft + # 2 - 2m/) | 



(21) 



graviton - scalar - scalar vertex 
5 



" f k lp k 2a C^-M 2 sV ^ 



(22) 



where S stands for the Goldstone (p^ of the gauge boson W. The first line is the contribution 
coming from the minimal energy-momentum tensor while the second is due to the improvement 
term. 
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graviton - scalar - fermion - fermion vertex 




Sipip 



Pl + C 



Sipip 



Pr) *T 



where the coefficients are defined as 

e m. 



(fi+ipip 

y <fr~ip-ip 



V2s w m w <«" 



—i 



V2s w m w 



V; 



Iplp 1 



\f2s w mw W 
graviton - gauge boson - fermion - fermion vertex 



C R - =i e ^V*- 




~{c^p l + c^p r ) C^f, 



with 



a 



w+i>ip 



-v. 



-VI 



C\ 



R 

W ± tpip 



0. 



(23) 



(24) 



(25) 



\/2sw ' ' y/'ZSw 
The tensor structures C, D and E which appear in the Feynman rules defined above are given by 

Cfivpa = Vpp Vvcr + Vpcr T]i/p 1]/iu T)pa j 



(26) 



Df Jll/pa (k 1 ,k2) = Vpu k la k 2p 



if°k!{k p 2 + 7]^ p k la k 2u - VP* h n k 2 „ + (fJ, o v) 



Ep l/ pa(ki,k 2 ) = i)^ (kipki a + k 2p k 2a + ki p k 2u ) - 



Vua ki fl k lp + i) vp k 2p k 2a + (fj, o u) 



(27) 
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B Scalar integrals 



In this Appendix we collect the definitions of the scalar integrals appearing in the computation of the 
matrix element. One-, two- and three- point functions are denoted respectively as Ao, £>o arm Co with 



A (m 2 ) = [d n l ' 
m z J 



I 2 — m? ' 



o 

1 



CoOi, (Pi ~ P2) 2 ,pl,ml,mj,ml) = [ d n l- 



iir 2 J (I 2 — m§)((I +pi) 2 — m\){(l + P2) 2 — m\) 



(28) 



C Form Factors 

Here we list the coefficients C\ which appear in the expansion of the form factors F±, F3, F4, Fj, F$, 
Fu defined in Eq. (fT71) . The remaining ones, as already mentioned, can be computed using the Ward 
identities in Eq. (fl~9]) . 

• Coefficients C\ entering in F\ 



C° = | 3 ( x d - x s) + 6(1 - xs) (xf - Qx w ) + 16 (xf - xw) (x/ + 2x w ) j , 
2xf — 3xs + 4:Xw + 3 



C\ 



6A 

2 



C\ = — g { x D (6xf + x s — 16x w — 2) + x D (xs (-10xj + 18x w + 3) + 12 (xfXw + x] 

8 A X JJ 

+ Xf — 2x 2 v ) — 32x w — 3) — x 2 D (x 2 s (20a;/ — A6xw — 6) + 2:5 (— 28x fXw — ixj (7xf + 4) 

+ hQxly + 54x w + 4) + 24 (x f - x w ) (x f + 2x w ) + 6x f + x 3 s - 26x w - l) 

+ xd (x s (8xf + 6xw + 3) — 2xs (2xfXw +%f (2x/ + 5) — Axly + xw) — 8 (x/ — xw) 

x (x f + 2x w ) - 3x% + I0x w ) - (1 - 2x s f \2x w (2x f + x s ) + (x s - 2x f f - 8x^j 

Cf = lx D (- (6x f +x s - I6x w - 2)) + x D (x s (-10s/ + 18x w + 3) + 12 (x f x w + x\ 

+ Xf — 2x w ) — 32xw — 3) + x D (x s {20xf — A&xw — 6) + xs (—28xfXw — &Xf (7xf + 4) 

+ bGxly + 54x w + 4) + 24 (x f - x w ) (x f + 2x w ) + 6x f + x% - 26x w - l) 

+ xd (x s (8xf + 6xw + 3) — 2xs (2xfXw + Xf (2xf + 5) — Ax^ + x\y) — 8 (xf — xw) 

x {x f + 2x w ) - 3x s + Wx w ) + (1 - 2x s ) 2 (2x w (2x f + x s ) + (x s - 2x f f - 8x 2 v ^j 

Cf = -^{x 2 D (3{8x f -3)x s - 2x f (lQxf + 32x w + 11) + 10x^ + 9) - 3x 2 s (8x f + 26x w + 3) 
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+ 4x s (-7x f x w +x f (13a;/ + 5) + 42xw + 3Ax w ) + 24 (6a;/ - 7) x 2 w + 92x f x w 

- 2x f (2x f + 1) (12a;/ - 1) + 9a;| - 96x^ - 68x w | , 

I = ^ {*/ (8x w (x 2 D +x s -2)- 2x 2 D + 34 - 2x s ~ 36a^ + l) 

+ 4x w (x 2 D (-3x s + Ax w + 5) + {x s - 3x w - 2) {3x s - 2x w - 1)) - 12xj (x s - 1) + 12a^ j , 

1 = ^ [ 8x w (-2x 2 D + 5a;/ + 4a; s - 2) - Ax 2 w (x 2 D (-2x f - 2x s + 3) + (x f + x s ) (6x f + x s ) 

- 5x f - 2x s ) - 2x w (-Ax) (x 2 D +x s -2)+3x f (x 2 D (2x s -3)- (x s - 2) x s ) +2(x 2 D - x 2 s ) 
x (x 2 D - x s ) + Ax)) +Xf(2x f -x s + 1) (-x 2 D + (x s -2x f ) 2 + Ax f ) - 16x^| , 

I = -^^2x 2 D (x s (x f (-Ax f + 32x w + 7)-8xw-2) + (2-6xf)x 2 s + 2(3-16xf)x 2 w 

+ Ax f (4x f - 9)x w + x f (2x f - 1) (8a;/ + 5) + 8x w + 1) + x% (Ax f - 1) + 163;^ (I0x f + 9x s 

- 9) - 4a;^ (34a; f x s + 2Ax) - 50a; / + 3x s (9x s - 16) + 24) + 8x w (- (x f + 8) x 2 s + (4 (x f 

- l)x f + 6)x s + 6x f - Ax) (x f + 2) +4x 3 s - 2) - (x s - 2x f ) 2 (2x f (x s - 1) - 8x)+3x 2 s 

- 4z s + 2)-64*kj, (29) 
Coefficients C\ entering in F 3 



6A 



l - x s (3x 2 D + 32a;/ - 20x w + 3) + 2x f (7x% + Ax w + 9) - Ax w (2x 2 D + Ax w + 3) 



8x) + 6x 2 s 



2"T j - x 2 D (xs (6x f + x s + I2x w + 4) - 4a;/ - 8x w - 3) + x% 

• ~ x D ) ( 



3 A (X 2 g ■( ; ; 

+ xs (Axs — 3) (2xf + xs + Axw) 



q 2 



-J^x 5 D (-2a;/ +x s + 2x w + 2) + x% (Ax w (2x f + x s + 1) - 2a;/xs + 8x) 



4A 2 (x D + x s ) 

- 2x f - x% + x s - lGxw) + x 3 D (2x w (2x f (x s + 3) + x% + x s + 2) + Ax)x s + 2x f x s + 12x) 

- 2x f - 8 (x s + 3) xIy - x% - 2x 2 s - hx s + 2) + x% (2x 2 s (-2a;/ + x w - 4) + x s (2x f (2x w 
+ 9) + 4a;/ - 2x w (Ax w + 9) + 5) - Ax f x w -2x } (2x f + 5) + x% + 8a;^ + 12a;w - l) 

+ x D (x% (-8x f + Ax w - 1) + 2x 2 s (x f (Ax w + 6) + Ax) - x w (8x w + 9) + 2) 

- 2x s (x / (10a; w - 3) + Wx) + (3- 20x w ) x w + l) + 2 (5 - 2x f ) x w - Ax f (x f + 2) + 2x% 
+ 8a;^) + x 6 D + x s (-8a;/ + Ax w + 1) - 2a; s (2a;/ (5x w + 2) + 10ar^ - 5x w (Ax w + 1)) 

+ x 2 s (8x f (x w + 2) + 8x) - 2x w (8x w + 9) - l) + 4 (x f - x w ) (x f + 2x w ) + 2x\ 
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+ x^-j 2 (z£> + i)x f - (x D - 1) (x D + x s ) I , 

Cf = - 4A2 ^ _ ^ I - x 5 D (-2x f +x s + 2x w + 2) + x 4 D (4x w (2x f + x s + 1) - 2a;jar s + 8a;/ 

- 2xf — Xg + xs — 16xw) + %d {—^xw (2a;/ (xs + 3) + x 2 s + x s + 2) - Ax 2 fXs — 2xfX S — 12a;/ 
+ 2x f + 8 (as + 3) + x| + 2a| + 5x s - 2) + x% (2x 2 s {-2x f + x w - 4) + x s {2x f (2x w 

+ 9) + Ax) - 2x w {Ax w + 9) + 5) - 2x s (2x w + 5) - Ax 2 f + x s + Ax w {2x w + 3) - l) 

+ x D {x% {8x f - Ax w + 1) - 2x% (x f {Ax w + 6) + Ax 2 - x w {8x w + 9) + 2) 

+ x s (x f (20x w - 6) + 20x 2 - AOxly + 6x w + 2) + Ax } (x w + 2) + Ax 2 - 2x% - 2x w {Ax w 

+ 5)) + x% + x 3 s {-8x f + Ax w + 1) - 2x s (2x f (5x w + 2) + lOx 2 - 5x w (Ax w + 1)) 

+ x 2 s (8x f (x w + 2) + 8x 2 - 2x w (8x w + 9) - l) +4 (a;/ - (a;/ + 2a; w ) + 2x1, j 

2g 2 f 

- X-rS 2 (xr> - l)xf + {x D + 1) (xd - xs) 



A 

Cl = -^^x 2 D (-12xs(2xf+x w )+Axf(3-8x w )-28x 2 f + 3x 2 s +Ax w (3x w + 5) + l) 

+ x% (8x f ~2)+x s (Ax f {Ax w - 7) + 44x/ + 4x w {9x w - 1) - l) + (22a;/ + 2) x 2 s 

- 8 (-i/iff (9x w + 2) + 3x 3 f + 2x 2 + §x 2 w (x w + 1)) + 10a;/ - 3x| - 4x w j , 

Cl = -^|-a;|,(4(a;/-10)a; w + 2a;/ + 15a;s(2a; V y + l)+44a;^-7)+a;|j(8a; w +4; 

+ x 2 s (-18a;/ + 28x w + 5) + x s (8x f (x w + 5) + 36a;/ + 26x w (2x w - 3) - 4) 

- A(x f (-Wxly + i^ + 5) + 6x 3 f + 9x 2 f + 2x w (<ox 2 w + x w - A)) + 3x| j 

+ X^y (x 2 D - 2x f - x s ) , 



-A 



Cl = ^^A0x 3 v (-x 2 D + 2xf+xs)-Ax 2 v (-x 2 D {8x f +Axs + 5) + 2x 4 D + A(x f + 2)x s 

+ A (3xj + x f - l) + 3x|) + 2x w (x 2 D - 2x f - x s ) (5x 2 D - 4 (x f + 3) x s + Ax f (x f + 1) + x% 
+ 6) + Ax 2 D x s -x%- 2x 2 D - 8x f Xg + 2Ax 2 x 2 s + 2Ax f x 2 s - 32x 3 x 5 - A8x 2 x s - 2Ax f x s 

+ 16xj + 32x 3 f + 2Ax 2 f + 8x f +x%- Ax% + 2x 2 s - 32x^ j 

" Z ^x 2 D (2x w - 1) - 2x w (2a;/ + x s ) + (xg - 2x f ) 2 + 4x/ j , 
-| - 2a; w (x^ + 2x f + x s - 2) + (-2a;/ + x s - 1) (a;^ + 2a;/ - x s ) + 8x 2 w ^ 
x |a;|, (4a;/ - 1) - 4x s (x f + x w ) + A ([x f - x w ) 2 + x w ^j + x|| , 

• Coefficients C\ entering in F4 



r 7 

° 3 " 4A2 



(30) 
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3^2 jz/ ( x d - 7«s + 10x w + 6) -x w (7x% - 19x s + 20x w + 12) + lO^j , 
3 g 2^2 (Jj _ x 2^ ^ x d (6xfX S - lOxf + 12x s x w - 2x% + x s - 20x w - 3) + 2x% 
(xs - 3) x s {2x f +x s + 4xw) | , 

A 3 ^ +Xs ^ j^D (^Z fas + 4a; w + 2) + 4x^ + x w (-9x s - 8x w + 4)) + x|, (x/ (2 (x s 

6) x w + x% - 2x s + 7) + 24 {xs + 6)-x w (4 (x s + 6) x w + 7x% - 10x s + 5)) 
2x 2 D {xf (2x s {-4xg + 5x w + 6) - 7x w - 4) + x^ (10x s - 7) + x w (5x s (3x s - 4x w - 4) 

14x w + 7)) + x D {-x f {x s {x s (12x s - I6x w - 15) + 22x w + 6) + 8x w + 3) 

2x) {x s (8x s - 11) - 4) + x w {x s {x s (24x s - 32x w - 29) + 44x w + 2) + 16x w + 7)) 

2xf,x w + x w {2x f (2 (x s - 4)x s + 1) + x s (6 (x s -2)x s + 7))+x f {x s (4x/ (x s - 4) 

2x 2 + 2x s - 3) + 2x f ) - 4 (2 (x s - 4) x s + 1) x^j , 

A 3 ^ I - x 4 D {x f (x 5 + 4:X W + 2) + Ax 2 + x w (-9x 5 - 8x w + 4)) + x 3 D {x f (2 {xg 

6) i ff +4- 2x s + 7) + 2x^. (x s + 6) - x w (4 (xs + 6) x w + 7x% - 10x s + 5)) 
2x 2 D {x f (2x s (4x 5 - 5x w - 6) + 7x w + 4) + x) (7 - 10x s ) + x w (5x 5 (-3x 5 + 4x w + 4) 
7 (2x w + 1))) + x D {-x f (x s (x s (12x s - 16xw - 15) + 22x w + 6) + 8x w + 3) 
2x^ (x s (8x s - 11) - 4) + x w {x s {x s (24x s - 32x w - 29) + 44x w + 2) + 16x w + 7)) 
2x^x w - x w {2x f (2 (x s - 4) x s + 1) + x s (6 (x s -2)x s + 7))+x f {x s (2x s {-2x f 

x s - 1) + 16x/ + 3) - 2x f ) + 4 (2 {x s - 4) x s + 1) x^j , 

^3 j^D ( x s (36a;/ - 6x w - 9) - 2 (2x/x w + 5x/ (8x/ + 3) + 6x^) + 6x| + 40x w + ll) 
4xjj {2x f + 1) - x% {lOxf + 108x w + 5) + 2x 5 (2 {-58x f x w + 2hx 2 +x f + 81x w ) + 77 x w 

I) - 4 (6(13- lhxf)x w - 59x f x w +Xf {5x f {6x f + 1) - 2) + 60x^, + 20x w ) +3x|| , 

^3 I - x 2 D {-3xg {8x f - 26x w + 3) + 50x/ (2x w + 1) + 24x^ + 6x| - 4x w {55x w + 4) 

II) + Ax 4 D (8x w + 1) - 20x^ (18a;/ + 13x s - 2) + 2 (2x s - 1) x w (50x/ + 7x s + 8) 
42x/x| - 132xjxs - 56x/xs + VlQx) + 156x^ + 40x/ - 3x| + 5x| + 2x s + 240x^| , 

^3 j^D { 2x w (8x/ - 4x 5 + 5) + Axf - 2x s - 56x w + 3) + 2x 2 D (x| (-6a:/ + 13a; w - 2) 

x s (-5 (4x/ + 5)x w + I2x 2 + 16a;/ + 52x^ + l) + 2 (5 (12a;/ + l)x^ + (5 - 6x f )x f x w 
x f {Axf {x s + 3) + 7) - 50x^) + x% + \0x w - l) - 2x% (8x/ + 5x w + 1) + x% {x f (28 
Ux w ) + 72xj + 4 (6 - llx w ) x w + 2) - 4xs (6x/ (5 — 7xw) + Xf (60x^/ — 2Qxw + 4) 
32x} +x w (10(1 - 5x w )x w + 3)) +8 (50x/x^ + 2 (l - 15x 2 f ) x 2 w - x f {2x f (5x/ +9) 
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+ 7)x w + x f (x f + 1) (2x f (5x/ + 4) + 1) - 20*4) + 4} , 

jx^ (4*/ (4x/ - 2x w + 1) + 2x s - 2x w - 3) - 2x 2 D (x 2 s (6x/ - 3x w - 2) 

+ x s (10a;/ (2a;/ - 4x w - 1) + 17xw + 1) + 6(4*/ - 3) x 2 w + 2x f (6a;/ + 25) x^/ - 4*/ (lOx 2 
+ x f -2)+x 3 s + 4x w - 16x w - l) + 8x^ (50a;/ + 37x s - 36) - \2x 2 w (8 (4a;/ - 3) x s 
+ 4(x f (5x f - 9) + 3) + 15x 2 s ) +2x w ((26x/ - 34) x| + (28a;/ (3a;/ - 2) + 22) x s + 44a;/ 

- 8X 2 - (5a;/ + 9) + 19x| - 8) - [x s - 2a;/) 2 (-20xj + (x s - 2) x<? + 2) - 160x^ J , 



Coefficients C\ entering in F-j 



... q 2x d 



(31) 
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6A 
2xi 



^3x 2 D + xs (4x/ — 4xw — 6) — 8 (xf — xw) (xj + 2xw) — 4a;/ + Axw + 3^ , 

- . 2 D — STI 2 ^ (~ 0*7 +a; s + 2xm/ - 1)) +x s (x s (-2a;/ + £5 - 4x w - 1) 
6 A (x D x s ) I 

+ 6(x/ + 2x w ))-3(x/ + 2x w ) j, 

7 = 4A2 ^ - I - x% (x f (-2x5 + 4x w + 2) + 4xj + 2x w (6x g - 4x w + 1) - l) 

+ x|j (3x| (4x/ — 4xvi/ — 1) + xg (34xw — 6 (AxfXw + 4x 2 + Xf — Sx 2 ^)) + 16 (x/ — xw) 

x (x/ + 2x w ) + 2x f - 18x w + l) + x% (x| (8x/ - 8x w - 2) + x| (-4a;/ (4xw + 1) - 16x/ 

+ Axw (8xw + 9) + 1) — 2x5 (3x/ + xw + 1) + 16 (x/ — xw) (xj + 2xw) + 6x/ — lOx^v + 1) 

+ x D (x s (x s (-24x/ + 8x w - 7) + 8 (x f - x w ) (x f + 2x w ) + 22x f + 6x| - 2x w + 2) 

- 2 (2x f + x w )) + x 5 D (-2a;/ - 8x w + 1) + (1 - 2x s ) 2 (-2x f + x s + 2x w ) (x s 

- 2 (x f + 2x w )) I - x^|2x/ (x D + 2x s - 1) - x D x s + x 2 D +x D - 2x 2 s + x 5 | , 

7 = 4A2 ^ _ xs - ) I - x% (x f (-2xg + 4x w + 2) + Ax} + 2x w (6x s - 4x w + 1) - l) 

+ x 3 D (3x| (— 4x/ + 4xw + 1) + x s (6 (AxfXw + 4x 2 + x/ — 8x w ) - 3Ax w ) — 16 (x/ — x w ) 

x (x/ + 2x w ) - 2x f + 18x w -l) + x 2 D (x% (8x/ - 8x w - 2) + x% (-4x/ (Ax w + 1) - 16x 2 

+ 4xh/ (8x w + 9) + 1) - 2x5 (3x/ + x w + 1) + 16 (x f - x w ) (x/ + 2x w ) + 6a;/ - 10x w + 1) 

+ x_d (2 (2x/ + x w ) - x s (xs (-24x/ + 8x w - 7) + 8 (x/ - Xvk) (a>/ + 2xw) + 22x/ + 6x| 

2x w + 2)) + x 6 D (2x f + 8x w - 1) + (1 - 2x 5 ) 2 (-2x/ + x s + 2x w ) (x s - 2 (x/ + 2x w )) J 

2<Z 2 f 

X— i (-2x/ + x s - 1) + x fl - (2x s - 1) (xs - 2x/) 
(- 



A 

q 2 x D 



6A 2 

- 4x/ + 3x| + 12x w + 1 



2x/ + xs — 4xvy — 1) |x 2 D (8x/ - 2) - 2x5 (2x/ + 6x w + 1) + 12 (x/ - x w f 
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Cf = ^^|-48a;^ + 4a;^.(-8a^ + 18a:/ + 7a;s + l) + (-2x/ + a;s-l)(a^-4(3a;/ + 2)a:fif 
+ 4{3x f (xf + l) + l) + 3xl)+2\x w (-8xf+4xs + 9)^+x^^-(2xf-x s + l) , 



C? = -^-^]2a;5,(-a:s(a; W '(8a;/+4!r w - + 23)-2)+a: w -(a;/(8x w --2) + 8a;J 



4A 2 

+ 2 (7 - 8xiv) x w + 9) + 2x 2 s x w -l)+x% (12x w - 1) - 2x| (4a;/ + 3a; w + 2) 

+ x s (4 (pxfXw + 6xf (xf + 1) + x w ) + A2xw + 2) — 4xs (2xj (xw + 6) + xj (6 — 4xn/ (a;w 

+ 1)) +8^ +3^ (2(5- 3x w )x w + 7)) +8 ((10a;/ + l)x%, - 2(x f + 1) (3a;/ - l)a;^ 



2g 4 a;D 
+ X — ; — 



x f (2x 2 f + x f + 2) x w + x f (x f + 1) {2x f (x f + 1) + 1) - 4a;^) + x% + 4x w j 
x l> + (#s - 2a;/) (-2a;/ + x s + 2x w ) + 4a;/ - 2x w j , 

4 r 

Cy = —^(-2x f + x s -4:X W -l)(2xf + xs-2x w -l)lx 2 D (4:Xf-l)-4:Xs(xf+x w ) 

+ 4((x/-x: w ) 2 +x: w ) +4}, (32) 



Coefficients C§ entering in Fg 



C& = ( x f - xw) (-2a;/ + x 5 - Ax w - 1) , 



Cl = 2WT~T 5T \ x d i Ax f - 5x s + 8x w + 5) +x s (5x s {-2x f + x s - Ax w - 1) 

6q A (x D x s ) I 

+ 12(x f + 2x w ))-6(x f + 2x w )^, 
c i = TT7~~ — ^^sa;i)(4a;/-a;s + 8a; w + 6)+a;c(9a;s(2a;/-a; s +4a;vF + 2)-16a;/-32a; w 



A 3 (x D + xs) 

- 7) - 2x 2 D (-2x 2 s (3 x f + Qx w + 2) + 9x f + 3x| + x s + 182:^ + 2) + x D (x s (x s (24a;/ 

- I2.T5 + 48^;^ + 17) - 10 (3a;/ + 6x w + 1)) + 4a;/ + 8x w + 1) + 2x 5 D 

- (2x 5 -l) 3 (xs-2(x/ + 2x: w )) J, 

C| = — ^ — <^ x^ (-Ax f + xs- Sx w -6) + x 3 D {9x s (2x / - x s + 4xw/ + 2) - 16a; / - 32xw 



A 3 (x_d - x s ) 

- 7) + 2a?|j (-2x% (3x f + 6a; w + 2) + 9x f + 3a;| + x s + 18a;iy + 2) + a; D (ar s (a?s (24a;/ 

- 12x s + A8x w + 17) - 10 (3a;/ + 6x w + 1)) + 4a;/ + 8x w + 1) + 2x 5 D 

+ (2.T S -l) 3 (a;s-2(a; / + 2 a ; H ,))|, 

Cf = g^j^D ((13 - 16a;/) x s + 32x/ (x/ + 2x w ) + 26a;/ - 34xvi/ - 13) + x% (2x f + 120x w + 19) 

- 2x s (-56a;/a; w + 2x f (x f + 2) + 150x^ + 86xw/ + 3) + 2 (-2a;/ (90x w + AAx w + l) 
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+ 60x 3 f - Ux 2 f + x w {3Qx w (4x w + 5) + 43) + l) - 15a;| j , 

Cf = ^j^w (-16^d + 90a;/ - 13x s + 29) + {-2x f +x s -l) (-13x 2 D - 4 (15a/ + 1) x s 

+ 60x f (x f + 1) + 15x% + 2) + 2\x w (-16a;/ + 8x s - 9) - 240^ j , 

Cf = ^fr{ 2a; l> C Ts ( 12 ( 2x f x w -Xf +xw) +7x w -4) - 6(4a;/ + 3)a;^ - 6a;/ (4a;/ + 5) x w 

+ I2x f (x f + 1) + x 2 s (3 - 6x w ) + 48a;^ - 5a; w + 2) + a;^ (4x w - 1) + 2a;| (20a;/ + 7a; w + 4) 

- 2x% (18a;/ + 2) + 60a;/ + (13 - 6x w ) x w + 2) + 4x s (-6a;/ (x w - 8) + 12a;/ (-3x^ 
+ x w + 1) + 40a;/ + x w (2 (a; w - 3) x w + 5)) + 4 (-2 (50a;/ + 13) x 3 w + 6 (2a;/ (5a;/ + 4) 

+ 1) x 2 w + (2x 2 (10a;/ + 9) - l) a;^ - 4a;/ (a;/ + 1) (5a;/ (x f + 1) + 1) + 40x^) - 5x| J , 

C| = ^^l2x 2 D (xs{-2x f {18x w + 5) + 15x w + 4:) + (6xf-3)x 2 s + 6(8xf-3)x 2 v 

- 3 (2a;/ (4a;/ - 9) + 5)x w - 2Ax) + 8x f - 2) + x 4 D (1 - 4a;/) - 2x% (2x f + 25x w + 4) 

+ 4a;s (72 (a;/ - l)a;^ - 3(2 (x f -3)x f + 5)x w +4x f - 2x f - 70x^) + x% (x f (8 - 36x w ) 

+ 90x w (2x w + 1) + 4) + 8 (5 (7 - 10a;/) x 3 w + 6 (x f - 1) (5a;/ - 3)x^ 

+ x f (2x f +3) (5a;/ - 3)x w + 2(2 - 5x f )x 3 f + 20x^) + 5x| + 20x w I , (33) 



• Coefficients entering in Fu 



2 f 

Oft = ^^|2a;/-3a;s + 28a;vK + 3 
^ = 3A (xj^- 4) ~ ^ + 6 ^ ~ XS ^ ( Xf + 2 ^ } ' 

°ii = „x2 V^, Ax i D {-2xf+x s + lQx w + 2)-x 3 D {x s {Qxf-AQx w -h) + 8{xf-x w ) 

4A^xd (xn + x S ) { 

x (a;/ + 2x w ) + x 2 s + 28x w + 2) - x 2 D (2x w (2x f (x s + 4) - 15x| + lla;s + 8) + 4x/x| 
+ Ax 2 x s - 8x/x,s + 16a;/ + 6a;/ - 8 (x s + 4) x^ + x% + l) + x D (x| (24a;/ - 2x w + 6) 

— 2x s (\QxfX W + Xf (10a;/ + 11) — 20x^ + 4x w + l) + 4 (xfx w + x 2 + Xf — 2x 2 v ) 

- 7x| - 2x w ) +x 5 D -{l- 2x s f (2x w (2x f + x s ) + (x s - 2x f ) 2 - 8x^) J , 

2 ( 

C " = -^^^^[- x D(- 2x f + x s + 10x w + 2)-x 3 D (xs(6x f -A0x w -5) 

+ 8 (x s - x w ) (x f + 2x w ) +x 2 s + 28x w + 2) + x 2 D (2x w (2x f (x s + 4) - 15x 2 s + llxs + 8) 
+ 4x f x 2 s + 4x 2 f x s - 8x f xs + 16a;/ + 6a;/ - 8 (x s + 4) x 2 w + x 3 s + l) + x D (x 2 s (24a; / - 2x w 
+ 6) - 2x s (l0x f x w + x f (10a;/ + 11) - 20x^ + Ax w + l) + 4 (xfX W + x 2 + x f - 2x 2 v ) 



17 



- 7x% - 2x w ) + x 5 D + (1 - 2x s f (2x w (2x f + x s ) + {x s - 2x f f - 8x^) 
C " = ^^^^ x D( 8x f- 2 )- 2x s(^x f + 15x w + l) + 12(x f ~4:X W )(x f -xw) 

- Ax } + 3x| + 2,Qx w + 1 
Cfi = ^rr-i -2x 2 D (8x w + l)-2x s (6x f -4x w + l) + 12(4xfX W +x 2 f +x f -5x 2 v ) 



3A 2 

+ 3x| + 8x w + 1 } , 



4 r 

Cfi = ^^l-12x 2 v (-2x 2 D + 6xf+xs + l)-2x w (x 2 D (4xf-2xs + 3) + (4xf-2)x s 
~ 4x f (3x f + 2)+ x 2 s ) + (-2x f +x s -l) (x 2 D - (x s - 2x f f -4x f ^j + iOx^ I , 



G\ Y = 1^^4x w (x 2 D (4xf-l) + 4xf(3x f +xs)-8xf + xs(3xs-4)+2)-(2xf+x s -l) 

x (x 2 D (4x f -l) + (x s -2x f ) 2 ^ -mx 2 w {2x s +x s -l) + 32a&, j . (34) 
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